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Abstract 



o 

(N 

$-H , This work is motivated by mathematical questions arising in differential equation models 

for autocatalytic reactions. In particular, this paper answers an open question posed by Guck- 
■<C ! enheimer and Scheper [SIAM J. Appl. Dyn. Syst. 10-1 (2011), pp. 92-128] and provides a 

' more general theoretical approach to parts of the work by Gucwa and Szmolyan [Discr. Cont. 

Dyn. Sys.-S. 2-4 (2009), pp. 783-806]. We extend the local theory of singularities in fast-slow 
polynomial vector fields to classes of unbounded manifolds which lose normal hyperbolicity due 
to an alignment of the tangent and normal bundles. A projective transformation is used to 
localize the unbounded problem. Then the blow-up method is employed to characterize the 
' loss of normal hyperbolicity for the transformed slow manifolds. Our analysis yields a rigorous 

■ scaling law for all unbounded manifolds which exhibit a power-law decay for the alignment with 
a fast subsystem domain. Furthermore, the proof also provides a technical extension of the 
blow-up method itself by augmenting the analysis with an optimality criterion for the blow-up 
exponents. 

(N : 

' Keywords: normal hyperbolicity, geometric singular perturbation theory, unbounded slow 

■ manifolds, blow-up method. 

a^ ■ 
p 

O ■ 1 Introduction 
(N , 

The motivation of this work are several models of singularly perturbed differential equations arising 
in applications. In [20j Merkin et al. propose to study a prototypical autocatalytic system given 
^ I by four reactions 

P-^A, A-^B, A + 2B-^3B, B-^C (1) 

where A, B are the two main reactants, P is a constant 'pool'-chemical and C is the product. Then 
it can be shown, using standard mass-action kinetics and non-dimensionalization [23], that ([I]) leads 
to a two-dimensional system of ordinary differential equations (ODEs) given by 

ex = yx^+y-x, 
y = fi-yx^-y, 

where ' denotes derivative with respect to time, x, y are dimensionless concentrations associated 
to B,A respectively and e,;U are parameters; we note that the parameter e is defined by the ratio 
of reaction rates for A ^ B and B ^ C \23\ p. 6192]. Obviously it is natural to assume that the 
concentrations are non-negative x,y £ Mq" := {w G R : w > 0}. Furthermore, note carefully that 
the nonlinear term arises due to the autocatalytic reaction part A + 2B — t- 3B. It has been proven 
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in [9] that the 2D-autocatalator ([2]) can exhibit an attracting relaxation-oscillation periodic orbit 
for certain ranges of the parameters; see also Figure [D^a). In [23] Petrov et al. generalized ([T]) by 
including a further reactant 

A, P + C ^ A + C, A^ B, B ^C, A + 2B^W, C ^ D. (3) 

As before, it is straightforward |23] to derive from ([3|) the ODEs 

ex = yx^ + y — X, 

y = fj.{K + z) - yx^ - y, (4) 
i = X — z, 

where (//,K,e) are parameters and {x,y,z) G (IRq^)^ are the concentrations. Numerical studies [HI 
[211 [22] have shown that periodic, mixed-mode and chaotic oscillations exist for the 3D autocatalator 




Figure 1: (a) Numerical illustration for the 2D autocatalator [2] with ^ = 1.1 and e = 0.01. The 
thin line segments show the vector field. The attracting periodic relaxation-type periodic orbit 
(thick black curve) is also shown, (b) The critical manifold Cq (gray) starts to align with the fast 
subsystem domains {y = const.} (dashed) as x — oo. 

For both models ([2|), Q it is often assumed that the ratio of time scales e > is a sufficiently 
small parameter; we shall also denote this assumption by < e ^ 1. In this case it follows that both 
autocatalator models are fast-slow (or singularly perturbed) ODEs. We use the notation p to denote 
a point in the phase space of concentrations for the autocatalator models i.e. p = {x, y) S (I^o^)^ or 
p = {x,y,z) G (K(j^)^. The x-nullcline, or critical manifold, for both models is given by 

= {p : yx^ + y - X = = \ p : y 



1 + x^ 

Note that Co is an unbounded smooth manifold. The angle between the tangent spaces TpCo and the 
hyperplanes {p '■ y = const.} decays to zero as x — )• +cxd; see Figure [l][|b). This alignment, which is 
expected to imply a loss of normal hyperbolicity in the system, already causes substantial difficulties 
in the rigorous analysis of the dynamics of the (2D) autocatalator model [9j. The global return 
mechanism induced by the unbounded part of Co also plays a key role for the complex oscillatory 
patterns of the 3D autocatalator [8, p. 98] which have been observed in numerical simulations. 

In a completely different context a model with similar properties to the autocatalator was 
proposed by Rankin et al. [24j as a caricature system for effects in aircraft ground dynamics 

e± = y + (x-/x)exp(f), 
ij = V — X, 
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where {fi,K,iy,e) are parameters, < e ^ 1 and Cq = {y = {fi — x)exp(x/K)} also aligns with 
{y = 0} as X — —oo if k > or as x — +oo if k < 0. 

Motivated by these examples, it is desirable to build a general theory of fast-slow systems with 
unbounded critical manifolds. The main contributions of this work are: 

• We study a general class of critical manifolds which may have an arbitrary power-law decay 
for the alignment with a fast subsystem domain. This includes the autocatalytic critical 
manifolds as special cases and answers open questions arising from various numerical studies. 

• Using the blow-up method we give a rigorous proof when normal hyperbolicity for a pertur- 
bation of the critical manifold fails. The relevant scaling law turns out to be given by 

(x, y) = (0(e-i/("+i)), 0(e^/(^+i))) , as e ^ 

where s is a power-law exponent computable from the critical manifold. This scaling has 
immediate consequences for the asymptotics of oscillatory patterns. 

• On a technical level we contribute to a further development of the blow-up method by adapting 
it to unbounded problems using projective transformations. Furthermore, we provide a novel 
argument for an 'optimality-criterion' of blow-up coefficients which have to be chosen in the 
analysis. 

The paper is structured as follows: In Section [2] we introduce the notation and the required 
background from fast-slow systems. In Section [3] a formal asymptotic argument is given to illustrate 
the important scaling properties and the system is localized via a projective transformation. The 
local system is desingularized via blow-up in Section [H The charts and transition functions between 
the two relevant charts are calculated. The dynamics in each chart are analyzed in Section [5] and [6] 
respectively. The optimality criterion for the blow-up coefficients is proven in Section [71 The main 
scaling result is summarized in Section [8] and its implications are discussed. 

2 Background and Notation 

In this paper we restrict our attention to the analysis of family of planar fast-slow systems given 
by 

f = y' = eg{x,y), 

where (x, y) G M^, < e ^ 1 and / : — M, : — M are assumed to be sufficiently smooth. 

Remark: By restricting to planar systems the (3D) autocatalator does not immediately fit 
within the theory we develop. However, we expect that center manifold techniques, similar to the 
generalization for folded singularities from three to arbitrary finite-dimensional fast-slow systems 
[25], can be applied. 

We are going to choose particular forms of /, g below but introduce some general terminology 
beforehand; for more detailed reviews/introductions to fast-slow systems see [121 [131 [2]. The critical 
set of ([6]) is given by 

Co:={{x,y)GR^:fix,y) = 0}. 
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We are going to assume that Cq defines a smooth manifold. Co is called normally hyperbolic at 
a point (x*,y*) G Cq if dxf{x*,y*) 7^ 0. We will assume that Cq is normally hyperbolic at every 
point in R-^ but Cq will be unbounded. Let us point out that there has been a lot of work on the 
geometric theory of planar fast-slow systems recently by De Maesschalck and Dumortier, see e.g. 
[171 flSl IT9] . particularly in the context of local singularities, periodic orbits and Lienard systems 
focusing on dynamics in fixed bounded domains. 

Remark: The general definition of normal hyperbolicity of an invariant manifold \10\ O [29] 
requires the splitting of tangent and normal dynamics. Recall that a compact manifold A4 C M'^ 
is normally hyperbolic for a vector field F : — >■ if there exists a continuous splitting of the 
tangent bundle 

TmR" = iv" e TM e 

where the linearization DF expands A^" and contracts A^* more sharply than TA4. In R^, we 
obviously just need one of the normal bundles if dim(A^) = 1. It is important to point out that 
we do restrict ourselves to perturbations of normally hyperbolic manifolds in the class of fast-slow 
vector fields of the form ([6|) being interested in loss of normal hyperbolicity of invariant manifolds 
within this class of perturbed vector fields where the perturbation is given by eg{x,y). 

The implicit function theorem yields that Cq is locally a graph. In our case, the parametrization 
will be global so that 

Co = {{x,y) :x = h{y)}. 

for some smooth function /i : R — )• R. The singular limit e — t- of ([6|) defines a differential algebraic 
equation, also called the fast subsystem 

x' = f(.x,y), , , 

y' = 0. 

Changing to the slow time scale r = et in ([6|) and taking the singular limit e = yields the slow 
subsystem defined on Cq 

= = f{x,y), 
^ = y = 9ix,y). ^ ^ 

The slow subsystem can be written as y = g{h{y),y). If h is invertible, the critical manifold can 
also be defined via Cq = {{x,y) G R^ : y = h''^{x) =: c(x)}. Implicit differentiation with respect to 
r yields y = c'{x)x so that the slow subsystem is 

c'{x)x = g{x, c{x)). 

Fenichel's Theorem [6l [T2| implies that that any compact submanifold A^o C Cq persists as a slow 
manifold Ad^ for < e <C 1. Furthermore, Ai^ is locally invariant and diffeomorphic to A^o- The 
flow on A^e converges to the slow flow as e — )• 0. For any fixed small e > we can define Ai^ and 
then extend it under the flow of ([6|) but this extension may no longer normally hyperbolic as an 
invariant manifold for the full system ([6|). To understand this effect we propose to study the family 
of model systems given by 

f{x,y,e) = l-x''y for s G N (9) 

and g{x, y) = fi for some ^ 7^ 0. The choice for the fast vector field is motivated by the asymptotic 
expansion of Cq for the autocatalator models ©J© as x — )• 00. In fact, returning to the chemical 
interpretation the general model class corresponds to autocatalytic reaction steps of the form A + 
{s + l)B ^ {s + 2)B. 



4 



The smooth unbounded critical manifold for Q is given by 



Co:=|(x,y)GM2.y=^=c(x)|. (10) 

As —7- oo the tangent space T(^x,y)Co starts to align with the x-axis since c{x) — t- and c'(x) — t- 
as |x| — )• oo. By reflection symmetry we restrict our focus to the positive part Cq H (M^)^ from now 
on i.e. the intersection with the non-negative quadrant will be understood. Note that the model 
class Q covers already a wide variety of power-law decay rates which seems like a good starting 
point for a general theory. Since we are only interested in the unbounded part of Cq we shall restrict 
the dynamics to 

n'' := {{x,y) eR^ :x>a} 

for a suitable fixed a > with a = 0(1) as e — )• 0; see also Figure [2]^a). Since dxf{x*,y*) = 
—s{x*y~^y* 7^ for {x*,y*) G Cq it follows that any compact submanifold A^o C Cq is normally 
hyperbolic and perturbs, by Fenichel's Theorem, to a normally hyperbolic slow manifold for e > 
sufficiently small. However, for fixed small e > we do not know yet when the extension of Ai^ 
under the flow starts to deviate from the critical manifold approximation. 

3 Asymptotic Expansion and Projective Transformation 

A good intuitive understanding of the problem is gained by considering one of the simplest possible 
systems in the model class Q given by s = 1 and ^ = 1 so that 

ex = 1 - yx, 
y = 1. 

From the viewpoint of formal asymptotics we can simply calculate a slow manifold by making 
the ansatz 

x{y) = XQ{y) + exi(y) + e^X2{y) H . (12) 

Inserting (|12|) into pi|) and collecting terms of different orders in e yields 

0{l): = l-yxQ{y), 

0{e'): ^ = -yxu{yl ^'^^ 
for A; G N. A direct calculation using (I13p gives the formal asymptotics of the slow manifold 

x{y) = \- f;(2fc - l)-i^ + 0(6^+i) as e ^ 0. (14) 
^ fc=i ^ 

The series ([H]) fails to be asymptotic when 1/y ~ iy2k+i as e — >• or y = 0{y/e). Since the 
critical manifold of (jlip is given as the graph of y = 1/x we expect that normal hyperbolicity of the 
associated slow manifold breaks down for x = 0(e~^/^). The example (jllh suggests that, depending 
on the asymptotics of the critical manifold as x — t- oo, there is an e-dependent scaling law that 
governs within which regime we can use the critical manifold approximation provided by Fenichel's 
Theorem. To give a rigorous proof of the asymptotic scaling we would like to desingularize the 
problem. In ^ the autocatalator is considered and phase space is re-scaled by an e-dependent 
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transformation; see also further remarks in Section [HI Here we use a different approach by first 
completely localizing the problem. Consider a projective coordinate transformation 

^:^1^(M2_^1) p{x,y)= (^^,y^ (15) 

where a = 1 is chosen for computational convenience for the domain T-L" . All the following results 
are easily checked to be independent of a fixed chosen a > independent of e; see also Figure [2][^b). 
A direct calculation yields the following result: 

Lemma 3.1. The map (|15p yields a transformed vector field defined on {M? — T-O') and given by 

ev = {yv^-'-v^), 
y = II. 




Figure 2: (a) Critical manifolds (jlOp for s = 1 (solid gray) and s = 2 (dashed black), (b) Critical 
manifolds after the transformation (llSp . Recall that we always restrict the dynamics to the non- 
negative quadrant (M^)^. 

The vector field (116p can be desingularized, similar to the desingularization for folded singular- 
ities [27] ■ via multiplication by and a time re-scaling which leads to 

ev = vHy-vn, (^7) 
y = fiv . 

Note that the transformation reverses time on orbits for t; < if s is odd but that we always restrict 
to the non-negative quadrant which implies that we do not have to consider this issue here. After 
the transformation the original critical manifold Cq = {y = 1/x'^} is given by 

p{Co) = {{v,y)€R^-n^:y = v'}=: So. 

Observe that Sq is normally hyperbolic except at the degenerate singularity {v, y) = (0, 0). The key 
point of the projective transformation p is that we have localized the approach towards a singular 
point to calculate the breakdown of normal hyperbolicity; in fact, it is not the singularity we are 
interested in but the scaling of the associated slow manifold 5^ in a sufficiently small ball excluding 
{v,y) = (0,0). Note that the case s = 1 yields the classical transcritical structure for Sq with 
f{v,y) = v{y — v); see [151 [26] aiid references therein. To desingularize one can try to apply 
a blow-up transformation. The blow-up method [3] was introduced into fast-slow systems theory 
in the work of Dumortier and Roussarie [3] . For additional background on the blow-up method for 
fast-slow local singularities see [16] . 
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4 Desingularization via Blow-Up 



Re-writing (fTTI) on the fast time scale and augmenting the system by e' = leads to the vector field 

v' = v'^ (y — v^) , 

y' = efiv', (18) 
e' = 0, 

which we denote by X. Consider the manifold B := x [0,ro] for some ro > where C 
denotes the unit sphere. Denote the coordinates on B by {v,y,e,f) where {v,y,e) £ S^. A general 
weighted blow-up transformation $ : ;B — t- M'^ is then given by 

V = f""-i;, y = r'^yy, e = r'^'e, (19) 

where we choose the coefficients as 

{ay,ay,ae) = {l,s,s + 1). (20) 

The choice of {av,ay,ae) can be motivated by using the Newton polygon (Sj H]. Via $ a vector 
field X is induced on the space B by requiring ^^{X) = X where is the usual push- forward. 
Recall that we restrict to the manifold Cq in the positive quadrant. Let By := B D {y > 0} and 
Bi := B n {e > 0} and consider two charts 

Ki:By^ and K2 : Bi 

where (ui,ri,ei) and {v2,y2-,f2) denote the respective coordinates for By and B^; see also [HI Fig2] 
for the basic geometry of a blow-up in M^. Slightly tedious, but straightforward, calculations yield 
the following result: 

Lemma 4.1. The maps ki and K2 are given by 

vi = vy"^^^, ri = fy-^/*, ei = ey~^^^^^/^, 

The transitions functions K12 and K21 between charts are 

V2=vie^ 2/2 = ei r2=rie^'^ 

-l/s 1/s -{s+l)/s 

VI = V2y2 , n = r2?/2 , ei = 2/2 ■ 



Lemma 4.2. The desingularized vector field in ki is given by 

v[ = svfil - vf) - fieivl~^\ 

r[ = firieivf, (21) 
The desingularized vector field in k,2 is given by 

¥2 = M^^2- 
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Proof, (of Lemma 14. 2p The required calculations follow by direct differentiation, for example, one 
has 

r[ = (fi / s)rl'^'^ eiv'l and v' = r[vi + riv[ 
which imply upon algebraic manipulation that 



s 

Similarly, one obtains the equation for ei given by 

/ + 1) s+l 2 s 

A division by the common factor rp~^ and a time re-scaling yield (j2ip . In K2 the transformation is 
just a re-scaling {v,y) = [€^^^^~^^'^V2,€^^^'^^^'^y2) and a time rescaling t — )• t/e gives (f22l) . □ 

5 Dynamics in the First Chart 

We can continue the critical manifold Sq = {{v,y) G (M^)"'" : y = v^} into the chart ki. Using the 
definition of the blow-up map and Lemma (|4.ip imply the relevant continuation result: 

Lemma 5.1. {ki o $-1)(5o) = {ivi,yi) G (M2)+ : vi = 1} =: 5i,a. 



Analyzing the flow (j2ip gives more information about Si^a- The vector field (j2ip has two 
invariant subspaces ei = and ri = 0. In the first subspace the flow is 

v[ = svl{l-vi), 

r[ = 0, ^^-^^ 

which has a line of equilibrium points Si^a which are attracting in the ui-direction. There is also a 
line of equilibria at ?;i = which we do not have to consider. In the second invariant subspace the 
flow is 

v[ = svl{l - vl) - ^lelv{+^ , , . 



Linearizing around (ui,ei) = (1,0) provides the important local dynamics: 

Lemma 5.2. The vector field (I24p has a center-stable equilibrium at {vi,ei) = (1,0) =: pi^a with 
eigenvalues and — and associated eigenvectors (— /i/s^,l)-^ and (1,0)^. 



Lemma 5.3. Locally near pa,i ihe equilibrium pi^a ho,s a one- dimensional center manifold 

Mi,a = {vi = l- ei^ + die? + 0(e?) := ni,,(ei) + 0(e?))} 
s 

where cu = -(1 -|- s -|- s^)/i^/(2s^). 
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Proof, (of Lemma 15. 3p Existence is guaranteed by the center manifold theorem [7] . Using a trans- 
lation Vi = vi — 1 niov6S pi.a to th.6 origin. A further linecir coordiiicite cli3,nge 

transforms (|24p into 

z[ = -fi{l + s)zl +0(3), 
4 = -s^ V(l + ^ + 2i^)^f + 2/uziZ2-(¥ + T)^i +0(3), 

where 0(3) = 0{zf,zfz2,ziZ2,Z2)- Making the ansatz Z2 = cuzf + 0{zf) and substituting into 
the invariance equation j7j for the center manifold at z = (0.0) implies the condition 

= (ens + ^^ + ll + ^^zl + 0{zf). 

Therefore, cn = —(1 + s + s'^) iJ? / {2s^) and transforming back to {vi, ei) yields the result. □ 

As before, other equilibria of (j24p will not be of relevance here. Lemma 15.21 and the center 
manifold theorem [7] imply the next result: 

Lemma 5.4. There exists a center-stable manifold Mi a for (|2ip atpi a containing Si and Ni^a- 
Furthermore, Mi^a is locally given as a graph of a map yi = hi^aifi,£i)- 

From the flow on the center manifold we see that there are two cases based upon the sign of 
fx. If /i > then trajectories in the center manifold flow away from the sphere while for < 0, 
trajectories in Aii^a How from the chart ki onto the sphere. We are only going to deal with the 
case fj, < from now on. The case > can be obtained from < by a time reversal of the the 
original problem (llSp . Figure [3l^a) provides a sketch of the results and the notation for the relevant 
dynamical objects. 

y2 /Af2,a 
n2,a{t) 

V2 



(b) 

Figure 3: (a) Sketch of the dynamics in chart ki. (b) Sketch of the dynamics in chart K2 
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6 Dynamics in the Second Chart 



Now the critical manifold in the chart ki has to be continued into the chart K2. This transition is 
studied using the continuation Ai2,a of 7Wi,a under the flow in K2; see also Figure [3jb). 

Lemma 6.1. The curve Ni^a C A^i,a transforms as 

Kl2{Afl,a) = |(w2,y2) G M2 : = 1 - + ^^SW + ^ (^2"^) 2/2 ^ +Oo| . 

We introduce the notation Af2,a ■= '^12(^/1, a) and note that A/'2,a can be extended under the flow 
()22p in chart K2 for initial conditions with 1^2 ~ 1 and arbitrarily large ^2 > 0. To emphasize that 
this extension is a trajectory we are going to denote it by n2^a{t) = {P'V2,ai^^i'^y2,a 
case distinction based upon the value of s will be necessary. For s = 2k with k €N the differential 
equation is 



v'2 = v^{y2-vf^ 
y'2 = M^'2^ 



(25) 



so that y2 < and V2 < for 1/2 < t^i'^ imply that nv2,a{t) — )■ for t ^ +cc The variational 
equation is defined via the time-dependent matrix 

vl{y2-vf) 

2nj,2,a(t)n„2,a(t) - (2/c + 2)n^.,^a{tf^^^ 1 
2A:/x n^a.aW^'^-^ 

which shows that the variational equation becomes asymptotically autonomous in forward time 
[25] and A{t) approaches a double zero eigenvalue when t — )• +00. Hence M.2,a is not everywhere 
normally hyperbolic in K2. A possible alternative argument for > 1, leading to the same con- 
clusion, involves desingularizing (I25p using division by V2 and observing that the resulting vector 
field is parallel on the 2/2-a'Xis with no flow in the y2-component. However, with this argument the 
case A: = is special but can be treated using standard results about the resulting Ricatti equation 
which also appears in the blow-up analysis of the fold point [14j . 

For s = 2A; + 1 and A; G No the vector field is 

v'2 = vl{y2-vf+^), 

y'2 = i^vf^'- ^ ^ 

Note that ()26p does not have the same nice monotonicity properties as before. To show the con- 
vergence towards ^2 = define 

y{v2.y2) ■■= + ^ §= Vf-'V'2 + ^,y2y'2 = -vf'^'^ < 

2k 2\fj,\ at 

showing that is a Lyapunov function [11] which implies the required convergence. Computing 
the variational equation yields 

A(t) - D, J ^2(2/2-^2' 

.- U{v2,y2) \ ,,„2fc+l 
V M^2 



"2,a(i) 

2n„2,a(t)%2,a(t) - {2k + 3)n^2,a(t)2'=+2 n^^^aitf 



^l{2k + l)n„,,„(t)2'= 
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which shows that the variational equation becomes asymptotically autonomous in forward time and 
A(t) approaches a double zero eigenvalue when t — t- +oo. As before, one could also use a suitable 
division by a power of V2 and observe that the resulting vector field is parallel on the y2-axis. 



Lemma 6.2. Ai2,a is not everywhere a normally hyperbolic manifold inside K2- 

Lemma 16.21 is just a local restatement of the alignment property of the original slow manifold 
Ce with the fast subsystem domains. This confirms the conjectured loss of normal hyperbolicity as 
X — )• +c». The previous results combine to give the following statement which is a partial version 
of our main result. 



Proposition 6.3. Consider the family of planar fast- slow systems 

x' = .27) 
y = efi. 

for s £ N, fi ^ and e > sufficiently small. Let A^o C Co be a compact submanifold of the critical 
manifold. Then Ai^ extends to a normally hyperbolic manifold of (|27p up to a domain of size 

(x, y) = (o(e-i/(^+i)), 0(e^/(^+i))) , as e ^ 0. 

Under the blow-up map (|19p - (|20p the manifold cannot be extended to a normally hyperbolic 
slow manifold to a subset of a larger domain. 

Proof. We work in projective coordinates via p. We can restrict to pieces of the critical manifold 
lying in the positive quadrant without loss of generality. Ad^ is obtained from A^o by Fenichel's 
Theorem since e > is sufficiently small. Using the blow-up transformation and Lemma 15. 4| it 
follows that Me extends to a normally hyperbolic manifold in the chart ki. By Lemma 16.21 the 
extension is not normally hyperbolic in the chart K2- By Lemma |4. II and Lemma [4.21 it follows that 
the neighborhood of {v, y) = (0, 0) scales as 

(.,y) = (ei/(^+i).2,6^/(^+^)y2). (28) 

By applying a blow-down transformation it follows that is normally hyperbolic up to a scaling 
region given by (I28p and not normally hyperbolic inside some larger domain. Via p^^ we get 
(x,7/) = (©(e"^/^*^^^), ©(e*/^*^^^)) as e —7- in original coordinates so that the result follows. □ 



7 Smaller Regions 

To prove the final result, we aim to strengthen Proposition 16.31 as it claims that the slow manifold 
extends and then normal hyperbolicity breaks down for the scaling (|28p using the blow-up (jl9p - 
()20p . A priori, we could have chosen different exponents for the blow-up which could have allowed 
us to extend the slow manifold even further. Here we show that the exponents in (j28p are indeed 
optimal in the sense that for any pair {ai,a2), > with < qi + a2 the manifold is not 
normally hyperbolic in the larger region 

(x, y) = (^0(e-(i+"i)/(*+i)), 0(e(^+"2)/(^+i))j ^ as e ^ 0. 
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Note that it suffices to assume that < ai + 02 ^ 1 so that the region is shghtly larger. Consider 
the modified blow-up 

V = fi+"iu, y = f'+°'^y, e = f"+^e. (29) 
Observe that in the chart K2 the blow-up (I29p reduces to the re-scaling 

which yields a strictly smaller scaling region than before. 
Lemma 7.1. In the chart ki the vector field is 

r> — _/f_^2-a2+s+sai s 

v[ = "^i^+"^V |-"^+^+^"^eir;f+^ + r|+"^+^t;^(r°^ - r^t^f), (30) 



Proposition 7.2. T/ie extension of the manifold 5^ in the chart ki for the blow-up (j29[l is noi 

normally hyperbolic as ri — )• 0. 

Proof. Suppose first that sai — 02 < 0. Since < ai + 02 *C 1 we can desingularize the vector 
field (j30p by 7.|^'^2+s+sai ]^ _ -)_ 5 -)_ sq;-|^ > g. Considering the invariant subspace defined by 
ei = of the desingularized vector field yields the equation 

which has a line of equilibria given by vi = r|"^ corresponding to the critical manifold Sq. 

Linearizing ()3ip around this line gives the variational equation 

which reduces to y = • y as ri — t- 0. Hence, in this case, the extension of is not normally 
hyperbolic in all of ki. Li the case sai — 02 > we can desingularize by where a* = 1 + 
s + min(Q;i + a2, sai — 0^2) in which case again a multiplicative factor containing ri appears in the 
variational equation. □ 

The main point of the previous proof is that we cannot remove the ri-dependence in the the 
chart Ki which means that there is no center manifold at the point pi^a extending up to the sphere 
5^, as previously for the blow-up with ai = = 02- 

8 The Main Result and Conclusions 

Finally, we can combine the previous result to obtain the main result about loss of normal hyper- 
bolicity and scaling. Recall again that we always restrict the dynamics to (M^)^ (iTi^ which is a 
subset in the non-negative quadrant bounded away from x = 0. 
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Theorem 8.1. Consider the family of planar fast-slow systems 



x' - 1 x'y, 
y = e/i. 

for s G N, /i 7^ and e > sufficiently small. Let d Cq he a compact submanifold of the critical 
manifold. Then Ai^ extends to a normally hyperbolic manifold of (I32p up to a domain of size 

{x, y) = (0(e-i/(^+i)), 0(e^/(^+i))) , as e ^ 0. (33) 

under the flow of (1321) . is not normally hyperbolic for any larger domain i.e. there does not 
exist a normally hyperbolic invariant manifold Ai^ in o- domain strictly larger than (j33p as e — t- 
such that is strictly contained in A4^. 

Proof. Apply Proposition 16.31 and Proposition 17. 2i □ 



Theorem 18.11 immediately applies to the 2D autocatalator. It is very interesting to point out 
the relation between the approach by Gucwa and Szmolyan [9j and the general Theorem 18.11 The 
scaling used in [9] to capture the large relaxation-type oscillation for the autocatalator is x = X/e. 
Theorem 18.11 states that the slow manifold loses normal hyperbolicity when x = 0{e^^^'^). This 
scaling certainly allows for the global crossings of the critical manifold when x = 0{e~^). In 
particular, the relaxation loop shown in Figure [T]^a) jumps near a fold point of the critical manifold 
'towards infinity' but it has a travel time so that it only turns around for x = 0{e~^). This yields 
one possible explanation why two consecutive blow-ups are needed in [9] as the scaling x = X/e 
compresses two important scaling regions for the global return simultaneously. 

For the model ([5]) by Rankin et al. the decay rate of the critical manifold is exponential and 
hence faster than any polynomial. In this context, Theorem 18.11 implies that a compact submanifold 
of the critical manifold will loose normal hyperbolicity for arbitrarily large choices of s S N in (|33p . 
Hence the scaling is expected to be x,y = 0(1), or at least close to this scaling for all practical 
purposes, so that normal hyperbolicity is lost in a very small region independent of e. 

As an important general conclusion from Theorem 18.11 we have proven a lower bound for large 
amplitude oscillations (LAOs) induced by a global return mechanism with the asymptotic dynamics 
given by p2p . For example, given a general autocatalytic reaction mechanism A -\- (s + 1)B — t- 
(s -|- 2)B then the lower bound is 0{e~^^^^^^^) for the fast component amplitude. Since s is usually 
known and e can be computed from the reaction rates the result has quite general applicability. 
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